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Measurement of J/ip Azimuthal Anisotropy in Au+Au Collisions at y/s NN = 200 GeV 
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The measurement of J/ip azimuthal anisotropy is presented as a function of transverse momentum 
for different centralities in Au+Au collisions at ^/s NN = 200 GeV. The measured J/tp elliptic flow 
is consistent with zero within errors for transverse momentum between 2 and 10 GeV/c. Our 



measurement suggests that J/ip with relatively large transverse momentum are not dominantly 
produced by coalescence from thermalized charm quarks, when comparing to model calculations. 



PACS numbers: 25.75.Cj, 12.38.Mh, 14.40. Pq 



Quantum chromodynamics (QCD) predicts a quark- 52 
gluon plasma (QGP) phase at extremely high tempera- 53 
ture and/or density, consisting of deconfincd quarks and 54 
gluons. Over the past twenty years, heavy quarkonia pro- 55 
duction in hot and dense nuclear matter has been a topic 5 e 
attracting growing interest. In rclativistic heavy-ion col- 57 
lisions the cc bound state is subject to dissociation due to 58 
the color screening and the change of binding potential 59 
in the deconfincd medium. As a consequence, the pro- 6 o 
duction of J/ip is expected to be suppressed, and such 6 i 
suppression has been proposed as a signature of QGP 62 
formation [1]. However, the J/ip suppression observed in 63 
experiments [2-6] can also be affected by additional cold 64 
[7, 8] and hot [9-14] nuclear effects. In particular the re- 65 
combination of J / ip from a thermalized charm quark and 66 
its antiquark has not been unambiguously established 67 
experimentally [11-14]. By measuring J/ip azimuthal 6S 
anisotropy, especially its second Fourier coefficient V2 (el- m 
liptic flow), one may disentangle J/ip from direct pQCD 70 
processes and from recombination. J/ip produced from n 
direct pQCD processes, which do not have initial collec- 72 
tive motion, should have little azimuthal preference. In 73 
non-central collisions, the produced J/ip will then gain 74 
limited azimuthal anisotropy from azimuthally different 75 
absorption due to the different path lengths in azimuth. 76 
On the other hand, J/ip produced from recombination of 77 
thermalized charm quarks will inherit the flow of charm 78 
quarks, exhibiting considerable flow. 79 

Many models that describe the experimental results ao 
of heavy-ion collisions depend on the assumption that ai 
light flavor quarks in the medium reach thermalization 82 
on a short timescale (~ 0.5 fm/c) [15, 16]. However, the 83 
dynamics through which the rapid thermalization hap- 84 
pens, are not very clear, and it has not been established as 
to what extent thermalization applies. The flow pattern as 
of heavy quarks provides a unique tool to test the ther- a 7 
malization. With much larger mass than that of light aa 
quarks, heavy quarks are more resistant to having their 89 
velocity changed, and are thus expected to thermalize 90 
much more slowly than light partons. If charm quarks 91 
are observed to have sizable collective motion, then light 92 
partons, which dominate the medium, should be fully 93 
thermalized. The charm quark flow can be measured 94 
through open [17] and closed charm particles. The J/ip 95 
is the most prominent for experiment among the latter, m 
However, because the J/ip production mechanism is not 97 
well understood, there is significant uncertainty associ- 98 
ated with this probe, since only J/ip from recombination 99 
of charm quarks inherit their flow. A detailed compari-100 
son between experimental measurements and models onioi 
J/ip i>2 vs. transverse momentum (pr) and centrality, inio2 



addition to nuclear modification factor, will shed light on 
the J/ip production mechanism and charm quark flow. 

This analysis benefits from a large amount of data 
taken during the RHIC [18] ^/snn = 200 GeV Au+Au 
run in the year 2010 by the new data acquisition sys- 
tem of STAR [19], capable of an event rate up to 1 kHz. 
In addition, the newly installed Time Of Flight (TOF) 
detector [20] allows STAR to improve electron identifica- 
tion, and background electrons from photon conversion 
are reduced by one order of magnitude due to less mate- 
rial around the center of the detector setup. 

The data presented consist of 360 million minimum 
bias (MB) events triggered by the coincidence of two 
Vertex Position Detectors [21], 270 million central events 
triggered by a large hit multiplicity in the TOF detec- 
tor [20] , and a set of high tower events triggered by sig- 
nals in the towers of Barrel Electromagnetic Calorimeter 
(BEMC) [22] exceeding certain thresholds (2.6, 3.5, 4.2, 
and 5.9 GeV). The high tower sample is equivalent to 
approximately 7 billion MB events for J /ip production in 
the high-pT region. In addition, in order to cope with 
the large data volume coming from collisions at high lu- 
minosity, a High Level Trigger (HLT) was implemented 
to reconstruct charged tracks online, select events with 
J/ip candidates and tag them for fast analysis. There are 
16 million J/ip enriched events selected by the HLT. 

The J/ip were reconstructed through the J/ip — > e + e~ 
channel, with a branching ratio of 5.9 %. The daughter- 
tracks of the J/ip were required to have more than 20 
hits in the Time Projection Chamber (TPC) [23], and 
a distance of closest approach less than 1 cm from the 
primary vertex. Low momentum electrons and positrons 
can be separated from hadrons by selecting on the in- 
verse velocity (0.97 < 1/(3 < 1.03), which is calculated 
from the time-of-flight measured by the TOF detector 
[20] and the path length measured by the TPC. At large 
momentum (p > 1.5 GeV/c), with the energy measured 
by towers from the BEMC [22], a cut of the momentum 
to energy ratio (0.3 < p/E < 1.5) was applied to select 
electrons and positrons. The electrons and positrons were 
then identified by their specific energy loss per unit track 
length ((dE/dx)) inside the TPC. More than 15 TPC 
hits were required to calculate (dE/dx). The (dE/dx) 
cut is asymmetric around the expected value for elec- 
tron, because the lower side is where the hadron (dE/dx) 
lies. It also varies according to whether the candidate 
track passes the 1/(3 and/or p/E cut to optimize effi- 
ciency and purity. The combination of cuts on 1//3, p/E 
and (dE/dx) enables electron/positron identification in 
a wide momentum range. Our measured J/ip particles 
cover the rapidity range — 1 < y < 1, favoring J/ip near 
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FIG. 1. (color online) Invariant mass spectrum of elec- 
tron/positron pairs for < pr < 2 GeV/c (a) and 6 < pr < 
10 GeV/c (b). The points are unlike-sign pairs with the J/ip 
signal. The blue solid line histogram shows the like-sign back- 155 
ground. 156 

157 
15S 

103 y = because of detection efficiency variation due to 159 

104 acceptance and decay kinematics. 160 

105 A total of just over 13000 J/ip were reconstructed in 1 " 
loe the entire pt range of 0-10 GeV/c, and the invariant massi 6 2 
107 spectra of e + e~ pairs with the J/ip signals for the lowest 
los (0-2 GeV/c) and highest (6-10 GeV/c) pr bins are shown 163 
109 in Fig. 1. 164 
no The following method has been used to calculate the 165 
in V2 of J/ip. Firstly, measurements of <p — 'F, ranging from 166 
n2 to 7r, were divided into 10 bins. Here <p is the az- 167 
u3 imuthal angle of J/ip, and <F is the azimuthal angle of 168 
u4 the event plane reconstructed from TPC tracks with the 169 
us azimuthally nonuniform detector efficiency corrected [24] . 170 
n6 Then two bins at supplementary angles were combined 171 
in into one. The J/ip yield within a combined <p — $ bin 172 
us was obtained by fitting the e + e~ pair invariant mass dis- 173 
u9 tribution with a Gaussian signal on top of a second order™ 

120 polynomial background. Then V2 was obtained by fit-175 

121 ting the J/ip yield vs. <p — 'F with a functional form of™ 

122 A{1 + 2v2Cos(2(</> — *)))• Finally, the observed was 177 

123 corrected for the event plane resolution [24]. ™ 

124 Three dominant sources of systematic error have been™ 

125 investigated for this measurement: assumptions in the u 2 lao 

126 calculation method, hadron contamination for the daugh-wi 

127 tcr e + e~ pairs, and the non-flow effect. The first source 1 ^ 

128 can be estimated from the difference in V2 calculated by 183 

129 methods with different assumptions. Two other meth-™ 

130 ods are used here. One is similar to the original method, 185 

131 except that the J/ip yield in each combined <p — 'F bin^ 

132 was not obtained from fitting, but from subtracting the^ 7 

133 like-sign background from unlike-sign distribution within 1 ^ 

134 the possible invariant mass range of J/ip- In the other 1 ^ 

135 method, the overall i>2 of both signal and background 190 

136 was measured first as a function of invariant mass, and" 1 

137 then it was fitted with an average of J/ip V2 and back-192 



ground V2 weighted by their respective yields vs. invari- 
ant mass [25]. The systematic error from hadron con- 
tamination can be estimated from the difference in cal- 
culated t>2 with different electron/positron identification 
cuts. While the original cuts aim for the best J/ip signal, 
a purer electron/positron sample can be obtained from 
a set of tighter cuts. The overall systematic uncertainty 
for the first two sources was estimated from the maxi- 
mum difference between the calculated V2 with the 3x2 
= 6 combinations of V2 methods and electron/positron 
identification cut sets mentioned above. Besides elliptic 
flow, there are also some other two- and many-particle 
correlations due to, for example, resonance decay and 
jet production. These non-flow correlations may influ- 
ence the reconstructed event plane and the measured V2- 
To estimate this non-flow influence on the V2 measure- 
ment, a method of scaling non-flow in p + p collisions 
to that in Au + Au collisions [26] was employed. This 
method assumes that 1) J/-0-hadron correlation in p + 
p collisions is entirely due to non-flow, and 2) the non- 
flow correlation to other particles per J/ip in Au + Au 
collisions is similar to that in p + p collisions. Under 
these assumptions, it can be deduced that the non-flow 
influence on measured J/ip V2 in Au+Au collisions is 
(V^ cos 2((pj/^, — <pi)) /Mv2- Here the sum is over hadrons 

i 

and the average is over J/ip in p + p collisions. M and 
are the multiplicity and average elliptic flow of hadrons 
in Au+Au collisions, respectively. Since the away side 
correlation may be greatly modified by the medium in 
heavy-ion collisions, this procedure gives an upper limit 
of the non-flow effect. Detector acceptance and efficiency 
variation with pt, centrality and rapidity may lead to a 
biased J/ip sample, which may induce some systematic 
effects when V2 also changes with these parameters. But 
these effects are estimated to be negligible compared to 
statistical errors. 

Figure 2 shows J / ip V2 as a function of transverse mo- 
mentum for different centralities, with the estimation of 
non-flow shown by the lines. Data from the central trig- 
ger, minimum bias trigger and high tower triggers are 
used for the 0-10 % most central bin, while only mini- 
mum bias and high tower triggered events are used for 
other centrality bins. Considering errors and the magni- 
tude of non-flow, J/ip V2 is consistent with for px > 2 
GeV/c for all measured centrality bins. Light particles 
usually have a larger V2 in the intermediate centrality 
than in the most central and peripheral collisions. This 
can be explained by a larger initial spatial eccentricity in 
the intermediate centrality, which is transferred into fi- 
nal state momentum anisotropy due to different pressure 
gradients in different directions, when there are sufficient 
interactions in the medium. However, no strong cen- 
trality dependence for J/ip V2 has been observed in our 
measurement. 

The top panel of Fig. 3 shows J/ip V2 for - 80 % 
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FIG. 2. (color online) V2 vs. pr for J/ip in different centrality 
bins. The brackets represent systematic errors estimated from 
differences between different methods and cuts. The pr bins 209 
for J/ip are 0-2, 2-5 and 5-10 GeV/c. The mean p T in each 210 
bin for the J/ip sample used for V2 calculation is drawn, but 2 " 
is shifted a little for some centralities so that all points can2i2 
be seen clearly. The curves show estimated non-flow innuence 2 i 3 
on measured V2 (see text). 
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FIG. 3. (color online) V2 vs. pr for J/ip in 0-80 % central 
events comparing with charged hadrons [27] and the <j> meson 
[28] (upper panel) and theoretical calculations [29-34] (lower 
panel). The brackets represent systematic errors estimated 
from differences between different methods and cuts. The pr 
bins for J/ip are 0-2, 2-4, 4-6 and 6-10 GeV/c, and the mean 
Pt in each bin for the J/ip sample used for V2 calculation is 
drawn. Estimated non-flow influence on measured V2 is shown 
as the shaded band. 



central collisions as a function of transverse momentum. 
For reference, two other sets of V2 measurements are also 
plotted, one is for charged hadrons (dominated by pions) 
[27] and the other is for the <p meson [28] which is heav- 
ier than the pion but not as heavy as the J/ip. Unlike V2 
of hadrons consisting of light quarks, J/ip V2 at pr > 2 
GeV/c is found to be consistent with zero within sta- 
tistical errors. However, the significant mass difference 
between J/tjj and light particles makes the direct com- 
parison of V2 vs. pt less conclusive. For example, for 
the same velocity at y = 0, the pr of J/ip at 3.0 GeV/c 
corresponds to px of pions (<p>) at 0.14 (1.0) GeV/c. 

In the bottom panel of Fig. 3, a comparison is made 
between the measured J/ip V2 and various theoretical cal- 
culations, and a quantitative level of difference is shown 
in Table I by % 2 /NDF and the p- value. «2 of J/ip pro- 
duced by initial pQCD processes is predicted to stay 
close to zero [29]. Although both normal suppression 
due to nuclear absorption and anomalous suppression in 
the hot medium due to color screening arc considered 
in the model, the azimuthally different suppression along 
the different path lengths in azimuth leads to a limited i>2 
beyond the sensitivity of the current measurement. On 
the contrary, if charm quarks get fully thermalized and 
J/ip are produced by coalescence from the thermalized 
flowing charm quarks at the freeze-out, the V2 of J/ip 
is predicted to reach almost the same maximum magni- 
tude as V2 of light flavor mesons, although at a larger pt 
(around 4 GeV/c) due to the significantly larger mass of 
J/ip [30]. This is 3a above the measurement for pr > 2 
GeV/c, leading to a large x 2 /NDF of 21.1/3 and a small 
p-valuc of 1.0 x 10 -4 , and is thus inconsistent with the 
data. Models that include J/ip from both initial produc- 
tion and coalescence production in the transport model 
[29, 31] predict a much smaller i> 2 [32, 33], and are con- 
sistent with our measurement. In these models, J/ip are 
formed continuously through the system evolution rather 
than at the freeze-out, so many J/ip could be formed from 
charm quarks whose V2 has still not fully developed. Fur- 
thermore, the initial production of J/ip with very limited 
V2 dominates at highp<r, thus the overall J/ip V2 does not 
rise rapidly as for light hadrons. This kind of model also 
describes the measured J/ip nuclear modification factor 
over a wide range of px and centrality [5] . The hydrody- 



TABLE I. Difference between model calculations and data. 
The p-value is the probability of observing a \ 2 that exceeds 
the current measured \ 2 by chance, even for a correct model. 

theoretical calculation 



X 2 /NDF p-value 



initially produced [29] 3.7 / 3 

coalescence from thermalized cc [30] 21.1 / 3 

initial + coalescence [32] 3.7 / 3 

initial + coalescence [33] 4.9 / 4 

hydrodynamic [34] 10.1 / 3 



2.9 x 10" 1 
1.0 x 10~ 4 
3.0 x 10" 1 
3.0 x 10" 1 
1.7 x 10 -2 
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256 
257 
258 
259 
260 
261 
262 
263 
264 
265 
266 



267 
26S 
269 
270 
271 



namic model, which assumes local thermal equilibrium, 272 
can be tuned to describe t>2 for light hadrons, but it pre- 273 
diets a J/ip V2 that rises strongly with p T in the region 274 
Pt < 4 GeV/c, and thus fails to describe the main fea- 275 
ture of the data [34]. For heavy particles such as J/ij), 
hydrodynamic predictions suffer from large uncertainties 27g 
related to viscous corrections (<$/) at freeze-out and ihe 279 
assumed freeze-out time or temperature. 230 

In summary, J ftp elliptic flow is presented as a func-J^ 
tion of transverse momentum for different centralities in 283 
y?SNN = 200 GeV Au+Au collisions. Unlike light flavor 
hadrons, J/ip vi at px > 2 GcV/c is consistent with zero285 
within statistical errors. Comparing to model calcula- 286 
tions, the measured J ftp V2 values disfavor the scenario 287 
that J/i[> with px > 2 GeV/c are produced dominantly' 
by coalescence from (anti-)charm quarks which are ther 
malizcd and flow with the medium. 
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